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Abstract 

We study limit distributions for random variables defined in terms of coefficients of a 
power series which is determined by a certain linear functional equation. Our technique 
combines the method of moments with the kernel method of algebraic combinatorics. As 
limiting distributions the area distributions of the Brownian excursion and meander oc- 
cur. As combinatorial applications we compute the area laws for discrete excursions and 
meanders with an arbitrary finite set of steps and the area distribution of column convex 
polyominoes. As a by-product of our approach we find the joint distribution of area and 
final altitude for meanders with an arbitrary step set, and for unconstrained Bernoulli 
walks (and hence for Brownian Motion) the joint distribution of signed areas and final 
altitude. We give these distributions in terms of their moments. 

Keywords: lattice path, catalytic variable, kernel method, moment method, Brownian 
motion area 



1 Introduction 

Banderier and Flajolet pQ give a nice treatment of the analytic combinatorics of directed 
lattice paths, i.e. walks on the integer line with jumps from a finite set. Their method to 
obtain the relevant generating functions is known as the kernel method [HE]. They address 
the asymptotic enumeration of unrestricted paths, bridges, meanders and excursions and 
further give limit laws for a couple of counting parameters on lattice paths, such as the final 
altitude or the number of contacts with the axis. In a later paper [2] the asymptotic study 
of the area under a path for excursions and meanders is initiated, by further using these 
techniques, however with the restriction to step sets containing a single negative step of unit 
length and the first moment only. This paper continues their work by completing the study of 
the full area distribution for arbitrary finite step sets. As limit distributions we find the area 
random variables of Brownian excursion and meander, which for some step sets is actually a 
simple consequence of weak convergence of the conditioned random walks to their Brownian 



motion counterparts, see the remarks below Theorem 2.1 In our combinatorial study we 
re-prove some of these results and provide the missing cases in the lattice path framework. 
Furthermore, following a suggestion in [23], we compute the joint distribution of the signed 
areas and the endpoint of a Brownian motion, which, to the best of our knowledge, are as yet 



unavailable in the literature. We obtain our results by studying a class of functional equations 



(3.3) which comprises those equations arising in lattice path and polygon counting [3]. This 
may be of use in different combinatorial contexts. 

Takacs \32\ IM] studies Bernoulli excursions and meanders to prove recursion formulas 
satisfied by the moments of the Brownian excursion area resp. meander area which have 
proved extremely useful in combinatorial probability. To that end he analyses a functional 
equation satisfied by the generating function of excursions E(z, q) (z marking length, q area) 
which reads E(z,q) = 1/(1 — zqE(zq,q). It allows to derive (univariate) generating functions 
for the excursion moments, see also |23| for a survey. This approach can be generalised to 
combinatorial classes counted by a "size" and an additional parameter, marked by z and 
q, whose generating functions satisfy equations of the form F(z,q) = G(z,q, F(zq,q)) [TOl 
I27j . This leads eventually to the very same recursions of the limiting moments and finally 
a Brownian excursion area limit distribution for that parameter. These "g-shifts" occur 
typically for cumulative counting parameters like construction costs of hash tables, internal 
path length in random trees or area of polynomioes, see e.g. j!3j . 

Such functional equations typically reflect a combinatorial decomposition of the class 
into smaller objects of the same class. For lattice paths with a more general set of steps 
than {±1} we are lacking a decomposition and functional equation as above. The kernel 
method has proved to be the weapon of choice for the enumeration by length in [lj. To 
that end, an additional "catalytic" counting parameter, the final altitude (marked by u) 
is introduced, which allows to turn a simple step- by-step construction of the lattice paths 
into a linear functional equation for the length-, area- and final altitude generating function 
F(z,q,u), however, with additional unknown functions in it. The shared feature with the 



above equation is the (/-shift, the occurrence of F(z,q,uq) in that equation, see (3.3) below. 
Very similar equations with (/-shift occur in different combinatorial contexts, e.g. column 
convex polyominoes counted by perimeter, area, and height of the rightmost column [3]. We 
prove a limit theorem for a class of functional equations which comprises the above mentioned 



lattice paths and also the column convex polygons in Section 11 



Remark. Upon the release of a preprint of this work, it has been brought to the author's 
attention that independently C. Banderier and B. Gittenberger have been continuing their 
work started in [2] intermittently over the past years, along similar lines. 



2 Statement of results 

We choose our notation close to the work in [1]. Given a finite subset 5cZwe consider 
walks on the integers with m steps from S, which are simply finite sequences (wq, w±, . . . , w m ) 
of integers with Wi — Wi-\ G S, i = 1, . . . , m. For the rest of the paper we fix wo = 0. The 
stress of this work is on walks (wq, . . . ,w m ), Wj G Z>o, in the non- negative half- line called 
meanders, and on meanders with w m = 0, called excursions. In order to make both problems 
non-void, we assume that there are positive integers c and d, such that 

— c = mine? and d = max S. 

The generating function of the step set is the Laurent polynomial 

S(u) = s- c u~ c + . . . + s d u d , 
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where Sj > 0, i = — c, . . . , d are non-negative weights (e.g. Sj G {0, 1}) and s_ c 7^ 7^ Sd- We 
shall occasionally speak of the mean and variance etc. of the step set, by which we mean the 
corresponding 5-valued random variable with probability generating function S(u)/S(l). 

For a step set S = {— c = o\ < 02 < • • • < 07 = d} we define its period to be p := 
gcd(<T2 — o"x, . . . , ai — 01). The step polynomial can hence be written as S(u) = u~ c H(u p ) with 
a polynomial H of degree (c + d)/p. We call 5 and the resulting walk aperiodic, if p = 1. 

The weight of a walk u> = (t«o, . . . , io m ) is the product 

Wt(u7) — S^j^ — togStoj — wi ' ' ' Sw m — w m -i 

of the weights of the steps of the walk. The "graph" {{i,Wi)\i = 0, ...,m} of a walk 
(wo, . . . ,Wm) can be viewed as a directed lattice path on Z x Z, whose steps are elements 
of {1} x S, meanders being directed paths never taking a step below the x-axis and excur- 
sions additionally ending on that axis, see Figure [TJ No confusion shall arise by identifying a 
meander with its graph. 




Figure 1: A meander and an excursion with steps from S = {—2, —1, 1}. 



For a meander w = (wq, . . . , w m ) with m steps we consider the functionals 

m 

l(w) = m, h(w) = w m , a(w) = Wj, 

i=Q 

which are the length, the terminal altitude and the area between the corresponding directed 
lattice path and the x-axis plus w m /2. Hence for excursions it is precisely the area. We refer 
to a in both cases as the area below a path. Our main result shows that this abuse of language 
is justified with view on asymptotics, as the expected area in the considered ensemble turns 
out to be at least of order m 3 / 2 while the final altitude w m is bounded by md. It proves useful 
for our approach to define the area random variables in terms of generating functions. Denote 
by F(z, q, u) be the generating function of meanders enumerated by the number of steps, area 
and final altitude, marked by z, q and u, respectively, i.e. the power series 

F{z,q,u) = J2^(w)z l{w) q a{w) u h{w \ 
w 

where w runs through the meanders with steps in S. Similarly, let Go(z,q) (= F(z, q, 0)) be 
the generating function of excursions. The area random variable X m for excursions and Z m 
for meanders is defined by 

\z m q l ]G (z,q) , \z m q l ] F(z,l,q) , 

F(X m = l)= \ , , P(Z m = l)= \ : , (2.1) {eq:XmZm} 

v ' [z m ]G (z, 1 ' y m > [z m ]F(z, 1,1) ' v ' x 4 f 
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where the square brackets denote coefficient extraction. Similarly, the random variable of 
terminal altitude H m on the set of meanders of length m is given by 

\z m u l ] F(z, u, 1) , x 

P <*" = '> = V]f(».U) ' (2 - 2) 

If, for example, all non-zero weights Sj are equal to 1, we consider the area and terminal 
altitude of a random meander (excursion) drawn uniformly from the set of all meanders 
(excursions) of length m. 



Brownian motion, excursion and meander 

The main results of this article is expressed in terms of integrals of processes related to 
Brownian motion. More precisely, let (B(t), t £ [0,1]) be a standard Brownian motion of 
duration 1, B + {t) = max{£>(t),0} and B~(t) = min{,6(i),0} its positive and negative part, 
respectively. We define the random variables 

A = f \B{t)\dt, A' = [ \B~(t)\dt, A + = [ B + (t)dt 
Jo Jo Jo 

of absolute, negative and positive area, respectively. The Brownian meander i3 me (i) (Brown- 
ian excursion B cx (t)) can be defined as the process (B(t), t £ [0, 1]) conditioned on B(t) > 0, 
t 6 [0,1] (on B(t) > and £>(1) = 0), cf. |llj . For a definition in terms of distribution 
functions we refer to |32l I34j . However, we only need their area random variables 



BMA= f B me {t)dt, B£A= f B ex (t)dt 
Jo Jo 



the (for our purposes) most convenient definition of which is in terms of their moments. 

Definition 2.1 (cf. [32]). The Brownian Excursion Area BE A is the random variable defined 
by the sequence of moments 

nr {btA j "^ r(3n/2-l/2)' ^ 
where the numbers K n are given by the recursion Kq = —1/2 and for n > 1 by 

3 — 4 n_1 

K n =—^K nr . 1 + ^2KiK n ^ l . (2.4) {eq:Kn} 

1=1 

Definition 2.2 (cf. [34]). The Brownian Meander Area BMA is the random variable defined 
by the sequence of moments 

±E(BMA") = ^ (1/ , 2)2 ""f N . (2.5) 
n! V ; Q r(3n/2 + 1/2)' v ' 

where the numbers Q n are given recursively by Qq = 1 and for n > 1 by 

3n — 2 n 
Qn = ^~ Qn~l+2Y,KiQ n -l- (2.6) {eq:Qn} 

1=1 



The numbers K n are given in (2.4). 
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Remark. BE A and BA4A are well defined by their moments, cf. |31j . 
With all the vocabulary at hand we state our results. 

Theorem 2.1. Let S be the step set of an aperiodic (see definition below) walk with step 
polynomial S(u) and drift 7 = S'(l)/5(1). Let furthermore r be the unique positive zero 
of S'(u) and (3 = \/2S{t)/ S"{t). The sequence of rescaled discrete meander area random 
variables of meanders with steps in S admits a limit law depending on the sign of the drift. 
More precisely we have for negative drift 7 < the weak limit 



pz m 

V2m 3 / 2 



B£A. 



(2.7) {eq:negdrif t} 



For zero drift 7 = the joint distribution of the rescaled area and terminal altitude random 
variable converges weakly to the joint distribution of (BA4A, B mc (l)), more precisely 



/3Z n 



y^m 3 / 2 ' v / 2m 1 /2 



(BMA, B me (l)). 



(2.8) {eq: zerodrif t} 



In both cases we have moment convergence. In the case of positive drift 7 > we find the 
expected area to be asymptotically equal to 7m 2 /2 while the standard deviation is of order 
o (m 2 ) . Hence the sequence is concentrated to the mean, 



Z n 



prob. 



7m 2 /2 



1. 



(2.9) {eq:posdrif t} 



The sequence of rescaled discrete excursion area random variables admits a weak limit law 
independent of 7, namely 



BSA. 



We also have moment convergence. 



Remark, i) Notice that in the case of zero drift r = 1 and /3/\/2 is the reciprocal of the 
standard deviation of the step set, in accordance with the rescaling of the random walks in 
the functional central limit theorem in [17j . Hence (2.8) follows from |17j by the continuity 
of the integral and evaluation operator / \-t /(l). 

ii) For 7 > the rescaled and "centered" discrete meander converges weakly to B [T7] and 
furthermore L B{i)dt is AA(0, l/3)-distributed. For the area we hence have 



7?n 2 /2 



where a 2 = [5"(1) + S'(l)] /5(1) 



am 



7 



3/2 



AA(0,l/3), 



is the variance of the step set and N(ijl,9 2 ) is the 



normal distribution with mean \i and variance 8 . 

Hi) Relations (2.7) and (2.10) do not seem to follow in this way: e.g. in [20] it is shown that 
a zero drift random walk of length n conditioned on having its first return to zero at time n 
converges weakly to the Brownian excursion. However, this conditioning does not imply that 
the walk is all positive or all negative prior to time n, since it may jump across without 
actually hitting (unless the step set is C { — 1, 0, 1}). Tying down a random walk as in [11J 
is not a continuous operation. Regarding negative drift meanders, we have only found weak 
convergence results for non-lattice step sets in the literature [12\ [21] . □ 



{eq: excursionla 



We formulate a more general result in Theorem 5.1, which Theorem 2.1 is a special case of. 
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Theorem 2.2. i) The moments of BE A allow the alternative description 



1 



E (B£A n 



-r(-i/2) 



2 ^n-l,l> 



n! v ' r(3n/2-l/2) 
where the numbers C n $ are defined by C n j = 0, if t < or n < 0, Co,o 

C n> t = C n: t-1 + (t + 2)C n _i j i + 2. 



(2.11) 



1 and efee by 



(2.12) 



Equivalently, C n ~ 



1,1 



8 n K n , with K n as in Definition 2.1 of BE A. 



ii) The joint distribution of the area BMA and the final altitude £> mo (l) of the Brownian 
Meander is uniquely determined by the sequence of moments 



M n>t = E [BMA n B 
where the numbers Q n j are given by 

Qn,t = < 

Qo,i 



l 



n!t!r(l/2)2-(™+*)/ 2 
r(3n/2 + t/2 + l/2) 



Qn,ti 



(2.13) 



Q„,t-2 + (i+l)Qn-l,t+l, if * > 1 
Qn-i,i + 2 • 8~ n C n -i i, if t = 0, 



(2.14) {eq:Qnt} 
1 and Q nj f = if t < or n < 0. In particular, we have 



with the initial values Qo,o 
Qn,o = Qn, cf. Definition 2.2, and Mot * s the tth moment of the Rayleigh distribution on 
[0, oo) given by the distribution function 1 — exp(— x 2 /2). 



Remark. The structure of the proof of the theorems |2.1| and 2.2 is as follows. From the 



functional equation (3.3) we derive moment convergence to limiting moments which satisfy 
recursion relations equivalent to those in Theorem 2.2 It is only in Section 10 where we make 
use of the fact that the Bernoulli excursions and meanders converge weakly to their Brownian 
siblings, in order to prove Theorem 2.2 Thus, assuming Theorem 2.2 is given, the limit laws 
in Theorem 2.1 can be viewed as a consequence of the functional equation only. 

The next theorem is about the joint distribution of signed areas and terminal altitude of a 
Brownian Motion. They follow from the convergence of the respective Bernoulli walks to their 
Brownian motion counterparts, cf. (9J [231 l32"j [34]. 

Theorem 2.3. The joint distribution of the signed areas A~ , A + and the terminal altitude 
B{1) is uniquely determined by its joint moments given by 



E 



{A+fiA-^Bilf 



L 



± 



(2.15) 



The numbers L,, t are given recursively by Lo,o,2t 
(M)^(0,0) by" 



r(3fc/2 + 3//2 + t/2 + 1) k ' l ' r 

1 and Lo,o,2H-i = 0, t > 0, and for 



r± 



k I 

3=1 3=1 



k,l—j,t t 2"'' u ^ fe '* ^ 2 ulc ' [) ^ l ' t ' 



(2.16) 



Similarly, the joint distribution of the absolute area A and B(l) is uniquely defined by the 
moment sequence 

raW2~( n+ *)/ 2 

E (A n B(lY) = w —j— 7 —^L n ,, (2.17) 



T(3n/2 + t/2 + l 
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where L n ^t+i = for n,t > 0, Lo 2t = 1 for t > and for n > 1 L n) 2t given by the recursion 

L n ,2t = 2 ^ KjL n ^j t 2t + Qn,2i- (2-18) 



XTie numbers K n are given in Definition\2.1\ and the numbers Q Uy t in Theorem 2.2 



3 The fundamental functional equation 

Our proofs are by the method of moments based on generating functions. Recall that the 
joint factorial moments of Z m and H m are given by 



E((Z m ) n (H m ) t ) = ' 



q=u=l 



[z m ]F(z,l,l) 

where (a) n = a{a— 1) • • • (a— n+1) denotes the falling factorial, similarly for X m . So, moments 
are expressed in terms of coefficients of univariate power series, which are in turn amenable to 
the process of singularity analysis. The required moment generating functions are "pumped" 
from a modified version of the fundamental functional equation of lattice path enumeration 

c-1 

F(z, q,u) = 1 + zS(uq)F(z, q, uq) - z ^ ri(uq)Gi{z, q). (3.1) 

where Gi(z,q) is the length and area generating function of meanders with terminal altitude 
i and r{(u) is the Laurent polynomial given by 

-(i+i)' 



n(u) =u l ^_ c u- c + ... + s_ (i+1) n-^J . (3.2) {eq:ri} 

This equation reflects the combinatorial decomposition "a meander is either i) empty or it is 
ii ) obtained by adding a step to a meander Hi ) without going below the x-axis" . The given 
functional equation is simply the g-shift of the fundamental functional equation for meanders 
in PQ, i.e. u is simply substituted by uq. The proof is similar to that in [2] where the authors 
study the actual area (without our notational abuse). 

Notation. For a commutative ring R and formal (commuting) variables oi, 02, ...,a n we 
denote by R[ai, 02, ... , a n ] the ring of polynomials, by R(a\,a2, ■ ■ ■ , a n ) the field of fractions 
thereof, i.e. rational functions, by i?[[ai, . . . , a n ]} the ring of formal power series. 

In what follows we study a slightly more general functional equation 

c-1 

F(z, q, u) = W(z, uq) + zS(z, uq)F(z, q, uq) — Tt(z, uq)Gi(z, q), (3.3) {eq:functionale> 

i=0 

where W(z,u), S(z,u) and the ri(z,u) are rational functions G M(z,u). Denote by Q(z,u) a 
fixed least common multiple of their denominators. We will have to take derivatives w.r.t. q 
and u of this equation, those w.r.t. q causing more trouble, since the multivariate chain rule 
[8J is involved. More precisely, for a function F(x\,X2) in two variables we need the formula 



\ „, , ^ fn\ n _ t ( d \ ( d 



\ n / \ 



dqj ' J^ \ t J \dx 2 J \dxi 



n—t / o \ t 



(F)(q,uq), 
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which is easily shown by induction. To make dealing with the derivatives of the functional 
equation a little less messy, we write for the derivatives of F(u, q) 



an+t 



(3.4) 

9=1 



and 

/ f> \ n 

(3.5) {eq:G_i} 

9=1 



:=G?\z) := (~) B G 4 (*,g) 



Notation: For bivariate functions T(z,y), we denote the partial derivative w.r.t. the second 
variable y by T'(z,y), T"(z,y), . . .,T^(z,y). 

An application of the above special case of the chain rule, setting q = 1 and regrouping 
terms eventually yields the following version of the nth ^-derivative of the functional equation 
evaluated at q = 1 : 

c-1 

{l-zS(z, u))F n >°(u) + n{z, ujG^ = zS(z, u)nF n - 1 ^ 1 (u) 

8=0 

+zS(z,u)f j (^jF^\n) 

+-ee (?) ( n t i y +t s {i) {^) F ^\u) 

c— 1 n / \ 

~zJ2Yl ( n i )u l r {l) (z,u)G < l n ~ l) +u n W (n) (z,u) =: RHS n {z,u). 
i=0 i=i ^ ' 



(3.6) {eq: derf cteq} 



If we multiply this equation by Q(z,u) the coefficients of F n '°(u) and the G,- on the left 
hand side are polynomials. It will turn out that for the asymptotic considerations only the 
terms in the top line of the previous equation are of interest. Further derivatives of equation 



(3.6) w.r.t. u yield an equation for F n,t (u). 

The following sections are dedicated to the singular behaviour of the functions G^\z), 
F n ' t {z, 1) and F n,t (z, ui(z)), where ui(z) is a branch of a certain algebraic equation. 



4 On the solution to the functional equation 



For a self-contained exposition, we sketch how equation (3.6) can be solved, assuming that 
the rhs is known. We then have c + 1 unknowns, and hence a seemingly underdetermined 
system. The kernel method [H El E] allows to solve this equation. 

We shall frequently deal with a c + 1-tuple (T(z,u), xo(z), . . . ,x c -i(z)) of formal power 
series T(z,u) G C[it][[z]] and xq(z), . . . , x c -i(z) £ C[[z]] defined by a linear functional equation 

c-1 

&(z, u)T(z, u) + ti( z , u)xi(z) = y(z, u), (4-1) {eq:generalf eq} 
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where &(z,u) and ti(z,u) are the polynomials 



$(z, it) = Q(z, it)(l - zS(z, it)), ti(z, u) = Q(z, u)n(z, u). 



In our framework y(z,u) € C[u][[z]] is a power series in z with its coefficients polynomials 
in it. The ti(z,u) are assumed to be linearly independent over C(z), the field of rational 
functions in z. The degree of <&(z, u) in u is assumed to be c + d with c, d > 0. We further 
assume that $(0, it) has degree c and hence there are c fractional power series u%(z), . . . , itc(^) 
and d fractional Laurent series vi(z), . . . ,Vd{z) with non-trivial principal part (each counted 
with multiplicities), such that §(z, Ui(z)) = = $(z, Vj(z)), see O Theorem 2] or [16]. 
We refer to the Ui{z) and Vj(z) as the small and large branches, respectively, and assume 
ui(z), . . . , u c (z) to be distinct. 

Definition 4.1. Let xq, . . . , x n be formal variables, and f{x) a formal power series. 
1. The divided differences are defined recursively as 



f [Xj] — / (Xj ) j / [%m , • • ■ , X 



/ \_ x n\ i ■ ■ ■ i %rik—i\ f [*^n2 ) • • ■ i ^"fc] 



Xni 



2. For the Vandermonde determinant in the variables xi, . . . ,x n we write 

A n (xi,...,x n ) = Yl (xi-Xj). 

l<i<j<n 

If /(x) is a formal power series in x, then f[xo, . . . formal power series in xo, ■ ■ ■ , x n , 

and furthermore, if Xq = . . . = x n then f[xo, . . . , x n ] = ^ f( n \xo). Furthermore, if / is a row 
vector of length n + 1, whose entries are formal power series, then the determinant of the 
n + 1 x n + 1 matrix with rows /(xj) factorises as 



det 



/ /(so) \ 
/(^i) 



det 



\ 



/KJ 

/[Xj , Xj-J 
/[Xi , • • • , Xj, ] 



V 



(Xj , • • • , Xj fc ) , 



/ 



where on the rhs the rows not involving Xj , . . . , xi k are unchanged. Notice that the determi- 
nant on the rhs is a symmetric function in Xj , . . . , Xi k [6]. 

Define the matrices 



M 



t (z,ui) 
t (z,u c ) 



t c -i(z,ui) 

tc-l(z,U c ) 



, J\f(z,u;y) 



y(z,u) t (z,u) 
y{z,ui) to(z,ui) 



y{z,u c ) t (z,u c ) 
We assume Ai to be invertible. By Ai we denote the matrix 

M = [tj(z, -)[iti, . . . ,Ui\) i=i,..., c 

j=0,...,c-l 



tc-l^jll) 
t c _i(z,ui) 

tc-l{z,U c ) 
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and by 



Gk(v) 



t (z,ui_) ... tjfc_i(z,ui) y(z, m) t fc+ i(z,ni) ... i c _i(z,ui) 
io^^c) ••• ife-l(^,«c) y(z,u c ) t k+1 (z,u c ) ... t c _i(z,u c ) 



and finally by .M(u;i) (resp. .M(it;i)) the matrix obtained from M. (resp. M) by substi- 
tuting by u. Notice that detAi and detAi(u;i) are the symmetric polynomials obtained 
by division of det M. (resp. Ai(u;i)) by the respective Vandermonde determinants. 



Lemma 4.2. The functional equation (4.1 ) has the unique power series solution with T(z, u) 
given by 



T(z,u) 



1 



$(z, u) 



det M(u;l) 
det M 



V 



<&(z, u) 

det M(z, u; y) 
®(z,u) detM 

and the power series x k {z) given 



-A det M(u;l) -A- 

(z,u)-^y(z, Ul (z)) H 



\ 



U — U-i 



1=1 



det M fj[ u i ~ u j 



(4.2) {eq: generalsol} 



Xk 



det g k (y) 
detM 



S-rZ 



^y{z,ui 



(j>k(ui, 



,Ul-l,U i+ l, 



.U c ) 



1=1 



detX 



n 



i 



(4.3) {eq:Gk} 



jJi u l ~ u 3 



where <fik( z ii ■ ■ • > z c— l) *s a symmetric polynomial in c — 1 variables. 

Proof. We drop the variable z and write y(ii) instead of y(z,u) etc.. Substituting u by 
«i, . . . ,u c eliminates T(z,u) from (4.1) and yields c equations for the Xi, namely the system 

M(x , . . . ,x c _i) T = (y(ui), . . .,y(u c )) T , 

By Cramer's rule and a Laplace expansion along the A;th column (k = 0, . . . , c — 1) we have 

det<7fc(2/(u)) _ <A , fc+ ;_ x det.Mfcj 



det M 



i=i 



(4.4) {eq:xk} 



where Aik,i is obtained from .M by deleting the Ith. row and the kth column. Since detM^j 
is an alternating polynomial in ux, . . . , ui-i, uz+i, . . . , u c , it can be factorised as 

det Mk,i = A c _i(ui, . . .,ui-i,ui+i, ... ,u c ) <fik(ui, . . . ,ui- 1: ui + i, ... ,u c ) 



with a symmetric polynomial (f>k [6], which implies (4.3). 
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As for the representation of T(z, u) we can write 



J2tk(u)x k = J2 t k(u)^^Y / y(ui)(-l) k+l - 1 detM k ,i 

k=0 fc=0 1=1 

E E «.(-)(-D w -' ** aiw = E ^i^^ 1 - 



(4.5) {alternative} 



fc=0 



This settles the first line of equation (4.2). The third line follows from a Laplace expansion 
of det Af (z , u; y) along the first column, and the second line from the □ 

Remark. The given representations of the Xk later allow to conveniently estimate the singular 



behaviours of the G 



(n) 



In Sections [6] and [7] we will need two further identities, whose proofs necessitate the following 
lemma. 

Lemma 4.3 (Grassmann-Plucker relations \22\ p. 34]). Let A £ c nxm annxm matrix with 
m < n. Let J, J' C {1, . . . , n} be index sets with \ J\ = \ J'\ = m and k € J\J'. Denote by Aj 
the submatrix consisting of the rows with indices in J and for t E J and an index I, Aj_ t+ i 
denotes the matrix obtained by replacing the row with index t by the one with index I in Aj. 
Then 

det ,4 j det ,4 j' = det Aj_ t+ i det Aj>-i +t . 

ieJ'\J 

□ 

Lemma 4.4. If we consider u,u±, . . . ,u c as formal variables, then we have the following 
identities 



d 

detMki tt- det J\(z, u; y(z, u)) 
' ou 



U=Ul 



d d 
det M t\ — det Gk - det Qk 7^ — det M , 

OU\ OU\ 



(4.6) {eq:Pluckerl} 



det M.(u;l) — — det Af(z, u; y(z, u)) 
ou 



U=U\ 



d d 
■ det M{z, u; y(z, u))-^ — det M. — detAi— — det At (z, u; y(z, u)). 



(4.7) {eq:Plucker2} 



Proof. On the lhs of (4.6 ) we apply two Laplace expansions to the (c+1) x (c+1) determinants. 



First expand along the leftmost column (containing y'(z,ui), y(z,Ui)) and then expand the 
resulting c x c determinants along the column containing the ife(iti). On the rhs expand each 
determinant along the row containing the tk(z,Ui) and y(z,Ui), respectively. The resulting 
representations are sums of terms of the form y(z, ui)tk(z, Uj) det Ax det A^ where Ax and 
Al are suitable submatrices of the (c + 1) x (c — 1) matrix 



A :-- 



t' (z,ni) ... t' k _ x (z,ux) t' k+1 (z,ui) ... t'^z^m) 
t (z,ui) ... t k _i(z,ui) £ fe+ i(z,iii) ... t c ^ 1 (z,m) 



t (z,u c 



t k -i(z,u c ) t k+ i(z,u c ) 



t c -i(z,u c ) 
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Let the rows be indexed with I = {— 1,1, ... ,c} from top to bottom, and for K C I denote by 
Ak the submatrix consisting of the rows indexed by K. We prove the identity by comparing the 
coefficients of y(z,ui)t k (z,Uj), y' (z,ui)t k (z,ui), y(z,ui)t' k (z,ui) and y'(z, ui)t' k (z, u\) j,l = 
1 . . . , c in both expansions. The coefficient of y(z, ui)t k (z, Uj), l,j £ I on the lhs is 



-l) l+ <» detA m>j} detAj\ { _ 1>1} 



and on the rhs 



-1 



(detA 



A{-i,i} det A i\{i,j} 



detA 



detA 



A{-i 



where e(j) G {0, 1} and we interpret the terms for j = —1 or Z = —1 accordingly. Let 
J = I\ {-1, 1}, J' = I\ {j, I}. If j, I > 2, then J\J> = {j, I}, J'\J = {— 1, 1}. If we apply 
the Grassmann-Pliicker relations with t = j G J \ J' the equality follows, observing that, for 
example, 



det 4r\{-i,z} det -4/\{ij} 



det Aj,_ { _ 1)+j det 



which is seen by elementary row operations. If j £ { — 1, 1} or / G { — 1, 1} one argues similarly 
(one of the summands on the rhs becomes zero) . 

To prove the identity (4.7), we expand the (c + 1) x (c + 1) determinants along the first 



column and compare the respective coefficients of y(z,Ui) and y'{z,ui) on both sides. For 
Y(z,ui) and Y'(z, U\) the equality is apparent, for % > 2 we apply the Grassmann-Pliicker 
relations to the (c + 2) x c matrix 



t' (z,ui) 
t (z,u) 
t (z,ui) 



i c _iO,it) 
t c -x(z,u\) 



t (z,U c ) ... t c -l(z,U c ) 

Here we we consider the rows indexed with the set / = { — 1, 0, 1, . . . , c} (from top to bottom), 
J = I \ {0, i}, J' = I \ {-1, 1} and t = -1. 

□ 



5 Evaluation at the small branches 

In this section we start the study of the singular behaviour of the functions F 11 ' 1 , which can be 
computed recursively with the method of the previous section from the functional equation 

c 

F(z, q, u) = W(z, uq) + zS(z, uq)F(z, q, uq) - r^z, uq)Gi(z, q). 

i=0 

The following assumptions imply that the singular behaviour of the F n,t and originates 
in that of a single small branch of the kernel equation. 
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5.1 Analytic assumptions 

Recall that $(z,u) = lc(z)(u — ui(z)) ■ • • (u — u c (z))(u — v\{z)) ■••(« — v c (z)) with lc(z) the 
coefficient of u c+d . We have to make some assumptions on the analytic behaviour of the 
branches Ui, Vj, or, equivalently, the singular points of the kernel equation §>(z, u) = 0. 

1. Small branches are distinct: <£(0,u) has degree c and the multiplicities of the small 
branches ui, . . . ,u c are one. 

2. Structural radius p and square root behaviour: There exists a point (p, r), < p, < r, 
with 1 - pS(p,r) = S'(p,r) = ( => $(p,r) = <&'(/>, r) = 0). Additionally we assume 
S"(p,r) > (i.e. &'(p,t) ^ 0). The equation z = l/S^-z, ii(z)) has hence two (real) 
solutions u (z) close to (z,u) = (p, r) with Puiseux expansions 



n ± (z)=r±/3 V / l-^/p + 0(l-z/ / 9) (*— >p), (3=^2^^. (5.1) {eq:gemil} 

We assume that u 1 * 1 can be continued analytically along (0, p) and that u~ has a finite 
limit at z = 0, while u + becomes infinite. Hence u~ is a small and u + a large branch, 
denoted by u\{z) and v\{z), respectively. 

3. Separation of small and large branches and uniqueness of p: We assume that z = p is 
the only point in {z : \z\ < p} for which u\{z) = v\{z). Furthermore 

V z G {w : \w\ < p} Vi G {2, . . . , c} V j G {2, . . . , d} : Ui(z) / vj(z). 

4. We assume that lc(z) ^ and Q(.z, u«(z)) ^ in {z : \z\ < p}. 

5. Zeroes of the determinant: Let g(z) be the greatest common divisor of the coefficients 
of t()(z, u), . . . , t c ^i(z, u) viewed as polynomials in u. We assume that 

Vz G {w : \w\ < p] : g(z)~ c detM ^ 0, 

or, equivalently, g(z)~ c det M. = 44> Uj(z) = Uj(z) for some i ^ j. 

6. We assume that 



lim |- (q(z,u)^(z,u) - ^ Q(z, m(*)) 



^0. 



Remark, i) The singular points of the equation <3?(z,u) = are to be sought for among the 
zeroes of \c{z) and of the resultant of &(z, u) and &(z, u) [16]. According to our assumptions, 
lc(0) = and p is a zero of the resultant. Furthermore, 4 is satisfied if the resultant of &(z, u) 
and Q(z,u) is non-zero in \z\ < p and lc(z) is non-zero in < \z\ < p. 

ii) Assumption 6 is in particular fulfilled if the degrees of Q(z, u)W(z, u) and the ti(z,u) 
as polynomials in u are at most c, since then at z = p the term in brackets is a polyno- 
mial of degree c with a simple zero u\{p) = r, and c — 1 not necessarily distinct zeroes 
u 2 (p), ■ ■ -,u c {p) ^ r. 

In the Meander case r is the unique positive zero of 5"(u), S"(u) > for u > 0, and 
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p = 1/S(t). Separation of the small and large branches follows from a domination prop- 
erty: For \z\ < p, i,j > 2 we have |"Uj(z)| < wi(|z|) and |^j(z)| > Ui(|z|). This and the 
uniqueness of p follow in turn from the additional assumption of aperiodicity of the step 
set, see pQ for details. Since lc(z) = s^z, all branches remain finite for z 7^ 0. Furthermore, 
Q(z,u) = u c = Q(z,u)W(z,u) and the degree of ti(z,u) equals i. Finally, g(z) = s- c z, and 
det M = (s_ c 2;) c A c (ui, . . . ,-u c ), i.e. detM = (s_ c z) c . 



We now state the general version of Theorem 2.1 where the conditions in items 1, 2 and 3 
generalise the drift being negative, zero and positive, respectively. 

Theorem 5.1. Let the sequences of (Z m ), (X m ) and (H m ) of discrete random variables be 
defined in terms of formal power series F(z,q,u) and Go(z,q) with non-negative coefficients 
as in (2.1) and (2.2), and let F(z, q, u) G C[q, u][[z]] and Gq(z, q), . . . , G c -i(z, q) G C[q] [[z\] 



be completely determined by afunctional equation (3.3), such that the above assumptions \5.1 
are fulfilled. 



1. If$(z,l) ^0 in {\z\ < p}, then (|2J]) holds. 

2. If T = 1 (i.e. §(p, 1) = 0) and $(z, 1) ^ in {\z\ <p}\ {p}, then (Q holds. 

3. If there is < z < p, ®{zq, 1) 
holds with 
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and ®(z, 1)^0 for {\z\ < z } \ {z }, then ( |279 ) 



S(z ,l) + z f z S(z,l) 



I2 = 2Q 



4- Equation (2.10) holds under either of the three conditions. 
Lemma 5.1. There is a > p such that the functions g\ n (z) are analytic in a slit disc 

D(p, a) = {z £ C | \z\ < a} \ [p, a) 
anrfF n,t (u) = F n ' t (z,u) is analytic in% = D(p,a)xC\{(z,Vj(z))\z G D(p,a), j = 1, . . . ,d}. 

Proof. By the above separation property we can choose a a such that the Ui(z) are separated 
from the Vj(z), i,j > 2, and z G D(p, a). By assumption 4. and general theory |16| a singular 
point zq G D(p,o~), zq ^ 0, for a small branch of the kernel equation can only be a finite 
branch point, i.e. two or more branches take the same finite value. Analytic continuation 
along a small circle with centre z$ permutes those branches coalescing there. 
Now we show that with y(z, u) analytic in % the function 



T(z,u) 



det Af(z, u; y) 
$(z,u) detM 



is analytic on %. To that end we first argue that T(z,u) assumes only finite values. Define 
the row vector 

R(z,w) = (y(z,w),t (z,w), . . .,t c -x{z,w)) , 

so the rows of J\f(z,u;y) are obtained by evaluating at w = u,u\, . . . ,u c . Let N(z,u;y) 
be the matrix obtained by replacement of the row R(z,Uj(z)) by the divided difference 
R(z, •) [u, u\, . . . , uj] , j = 1, . . . , c. Since R is analytic, these take a finite value at (z, u) = 
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(zq,uq), wether or not some of the branches coalesce in (zq,uq), by the continuity of the 
divided differences at multiple nodes. Then 



T(z,u) = ^Af(z,u;y) {eq:detNtilde} 

Hz)i\Ui( u - v j) detM 

since det AT(z, u; y) = A c+ i(u, m, . . . , u c ) det M{z, u; y). We hence have continuity at (zq, ^o)- 
By Hartogs' Theorem [29], analyticity of T(z, u) at (zq, uq) follows from the analyticity of 
the partial functions 

z i— )■ T(z, uq), u h->- T(zq,u). 

The latter is a meromorphic function in u, which has at most a removable singularity at no- 
For the analyticity of the former observe that the expression for T is a symmetric function in 
ui, . . . , u c . Analytic continuation of T(z, uq) along a small circle with centre zq permutes the 
branches Ui and takes T(z,uq) to itself. With Morera's theorem analyticity follows |16j . 

To prove the Lemma, one proceeds inductively. With y(z,u) = Q(z, u)W(z, u), we im- 
mediately obtain the assertions for the F 0,t (u), t > 0. For G,' (z) we can argue similarly 
with their respective determinantal representations. Now let the assertions be true for all 
(re — with n > 1. The determinantal representations of G k n \z) and F n '°(u), involve by 



(3.6) functions Gu\z) and F l,t {u), I < n — 1, I + t < n, which are analytic by the induction 



hypothesis. With the above reasoning the proof is complete. □ 

The qualitative singular behaviour at z = p is treated in the following lemma. 

Lemma 5.2. i) The functions G^\z), k = 0, . . . , c — 1, allow a Puiseux expansion in as- 
cending powers of yT — z/p. The same is true for F n,t {ui{z)) ifm(z) is analytic at z = p. 
ii) Assume that for an index set L C {2,...,c}, the branches ui(z), I G L form a cycle 
of length k at z = p, i.e. are conjugates of each other. Then F n,t (ui(z)), I £ L, allows a 
Puiseux expansion in ascending powers of (1 — zj p) l l 2p , if p is odd, and in ascending powers 
of (1 — z/p) l / p , if p is even. 

Proof. By the closure properties of algebraic functions and the algebraicity of the Uj, the 
functions in question are algebraic as rational functions in u±, . . . , u c - We can wr ite G^\z) = 
P(ui, . . . , u c ) and F n,t {u) = Q{u, u±, . . . , u c , where P and Q are symmetric in u\, . . . , u c . Now 
consider a small circle C with centre p. If we start at a zq G C and continue G^\z) resp. 
F n,t (ui) analytically describing C once counterclockwise, G^\zq) = P(u\(zq), . . . ,u c {zq)) is 
taken to P(vi(zq), ^2(^0), ••• , u c (zo)) and, if Ui is analytic at p, F n,t {u,j) is taken to 
Q(ui(zo) , v\(zq) , U2(zo) , . . . ,u c (zq)), since u\ is taken to v%, U{ remains fixed and the remain- 
ing branches are permuted. Describing C once more , we end up at Gj, n \z) resp. F n,t {ui) 
again. By Puiseux' Theorem [IB] , assertion i) follows. Let {w\, . . . , w p } = {ui, I £ L] in such 
a way, that describing C once counterclockwise takes Wj(z) into Wj+i(z) (indices mod p). 
Then describing C once takes F n,t {wj{zo)) to Q(vUj+i(zo), Vi(zq), uz(zo), ■ ■ ■ , u c (zo)). Notice 
that, for odd p, repeated analytic continuation of the pair (ui,wi) along C yields the cycle 
of length 2p 

(u^Wl) > (V!,W2) > {ui,W 3 ) > ... > (vx,W k -l) ► (ui,W k ) 

— > (vi,wi) — > (ui,w 2 ) — > ... — > (vi,w k ) — ► (ut,wi), 
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hence again by Puiseux' Theorem, the first assertion in ii) follows. In a similar fashion, 
for even p we have two disjoint cycles of length k, one containing all the (ui,v)2j) and one 
containing all the (ux, W2j+x)- This proves the second assertion of ii). □ 



(5.3) {eq:Sprime} 



By Tailor's formula the Puiseux expansion of S'(z,ui(z)) at z = p starts with 
S'(z, ui(z)) = S' (z, t - P^l-z/p + 0(1 - z/p)) 

= -S"(p,T)(3y^Z~/p~+0(l-z/p) 

Let similarly the small branches w\(z), . . . , w p {z) form a cycle at z = p, with a local expansion 

w j (z) = a + 1 ^- 1 \l-z/p) 1 /P + ... (*—►/>), (5.4) {eq:wj} 

where £ = exp(2iri/p) is a pth root of unity. Then we have 

$(p, a ) = $'(p,«) = •■■ = <S> {p ~ 1] (p,a) = 0, $ (p) (p,«) / 0, 
and the Puiseux expansion of &( l \z,Wj(z)) starts with 

$ W (z, wj{z)) = l _ $fr ) foa)7 P ~'£ (j '~ 1)Cp ~°(l - z/p) (p - l V p + ... (5.5) {eq:Phiwj} 

[pi). 

The functions F n,t (u\{z)), F n,t (l) and G , [™' > are easily seen to be rational functions in the 
Ui , and by the discussion above they admit a Laurent series expansion in powers of \J\ — z/p. 
We compute the leading singular term of these series, starting with F n,t {u\(z)). 

Lemma 5.3. The function F n,t (z, u) evaluated at the largest of the small branches u\(z) has 
the following leading singular behaviour at z = p : 

^(^M*)) = (1 _ ^/a+t/2+1/2 + ° ( (1 _ z/p )3n/2 + t/2 ) (* ^> P)> ^ 
where the numbers a n> t satisfy the recursion 

t 

2^' an ' t - 1 2 4 + 1 



a n ,t = tte ■ a n ,t-l + 7T " ■ i -, " fl n-l,f+2 (5.7) {eq:recursionan 



with 

£- (det M(z,u;Q-W))\ 

ao,0= hm — L /0 5.8 {eq:a00} 

z — >p 2Q{z,ui{z)) det M. 

In particular, this limit exists. If the small branches W\(z), . . . ,w p (z) form a cycle of length 
p at z = p with an expansion (5.4), then we have 

^(z, Wj (z)) = 0( (i _^ )3 1 n/2 _ 1/2+Vp ) (*-►,). (5.9) 

Ifui(z), i > 2, is regular at z = p, then 

F^(z,n t (z)) = o( (i _^ 1 )3n/2 _ 1/2 ) (*-►,). (5.10) 

Furthermore, for n > 1 £/ie singular behaviour of the functions G^ 1 ' (z) is 

4 n) (^)=0( (1 _ 2/p 1 )3w/2 - 1/2 ) (*->/>)• (5-11) {eq:singG k n} 
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Proof. If we consider the -F n,t organised in an array (n indexing rows, t columns), we can 
compute the entries inductively, row by row, where each skipping to the next row requires an 
application of the kernel method. We prove the lemma accordingly. 

First of all, ao,o is well defined by the continuity of the defining term at z = p, since a 
factor g(z) c A(ui, . . . , u c ) can be cancelled in the numerator and the denominator, and by 
Assumption 6 it is non-zero. 

Assertions for F°'°(ui) and F°''(uj), t > 0,i > 2: The assertion for (n,i) = (0,0) follows 
from the explicit form (TOJ of F°'°(u), where y(z,u) = Q(z,u)W(z,u). We have 



F°' u (ui(z)) 



l(detN(z,u;Q.W))\ u=ui 



ao,o 



-zQ(z,ui(z)S'(z,Ui(z)))detM - zip 



P), 



(5.12) 



where we used the fact that for a solution Uj of the kernel equation we have &'(z,Ui(z)) 
!))S'(z, Ui (z)). Fv 
1/S(p,r) and j3 



—zQ(z,Ui(z))S'(z,Ui(z)). Furthermore we used the expansion (5.3) and the defining rela- 

= y/2S(p,T)/S»( P ,T) 

(0) 



tions p 



Q(z,u)W(z,u) we see that G^' is finite at z 



Moreover, from (4.4) with y(z,u) = 

p. By the previous lemma it has at most 
a square root singularity at p. 

As for F°' t (uj(z)), j > 2,t > we look at a representation ( |5.2[ ) (and derivatives w.r.t. u 
thereof) to see that evaluation at u = uj(z) leads to a finite value at z = p for every t > 0. By 
the separation property and together with the previous lemma we obtain the desired bound. 
Induction step for u\(z) and branches ui(z) analytic at z = p: Let n > and assume 
that the estimate for and all F r,s (u\(z)) (resp. F r ' s (ui(z))) are established, where 

r < n — 1. We can compute the estimate for F n,t (ui(z)) for t > from the t + 1st derivative 
of (3.6) which can be written as (with the obvious adjustments in the case n = 0) 

(1 - zS{z, u))F n ' t+1 (u) - (t + l)zS'(z, u)F n '\u) = 

't + i 

2 

+ Rn,t+l(u) 



zS"(z, u)F n ^- 1 (u) + zS(z, u)nF n ^ t+2 (u) 



(5.13) {eq:Fnt} 



where R ny t+\{u) contains all the remaining derivatives of total order < n + t + 1 and of order 
< n in q. For u = Ui(z), we see by induction hypothesis that 



R 



n,t+l 



( Ul (z)) = (l-z/p) 



-(3n/2+t/2-l/2) 



and is hence of neglegible order, as z 

*)) 

F^\ Ul (z)) 



p. With the the above Tailor expansion (5.3) of 



S' (u\{z)) equation (5.13) reads on the level of leading order terms 

t a n ,t-i 



+ 



n 



2f3yfT^zJp~ (1 - Z /pfn/2+t/2 

S(p,r) a n -x,t+2 



+ 



t + 1 S"(p, T)(3yfT^zJp~ (1 - Z /p)3"/2+t/2 ^ (1 _ z//0 )3n/2+ t /2 



1 



P)- 



Recalling the definition of (5 this establishes the recursion equation. Notice that this settles 
in particular the assertions for F 0,t (ui(z). 

One argues similarly for the branches ui, where we first see that 



JM«K*)) = o ((i - z/pyw-w) 
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(with the G^ the dominating terms) and since S'(z,ui(z)) — > e / 0, we have the same 
estimate for F n '°(ui(z)). The estimate for F n,t (ui(z)) follows by induction on t. 

(n) 

Induction step n — 1 — )■ n for G k : Now let n > 1 and assume the assertions to be true 
for all (k,t) with k < n — 1. We apply now the kernel method to compute the , i.e. 
we substitute u\, ... ,u c into ( |3.6| ) (multiplied by Q(z,u)) and obtain c linear equations for 
G^ k , k = 0, ... , c — 1, the solution to which is given by (4.4), when y(z,u) is replaced by 
RHS n (z, u) = Q(z, u)RHS n (z, u). 

\k+c-l c 



G 



(n) 



-1 



n ) 



i=i 



^kjui, . . . ,Ul-l,Ul + l, . . . ,U C 

detM 



n 



By induction hypothesis, we have the estimates 

RHS n (z, Ui (z)) = 



(1 - z/p) 3n / 2 - 2+5 W 



(5.14) {eq:RHS} 



where 5(1) = 3/2, <5(i) = 0, if Ui(z) is analytic at z = p, and <5(i) = 1/p, if Uj(z) = Wj(z) 
as in ( |5.4[ ). Therefore, each summand is O ((1 — z/p)^ 3 ™/ 2-1 / 2 )) , since the products in the 
denominator is tends to a constant ^ for u\ and the branches ui analytic at z = p, while 
for a branch wj as in 

fliT4| ) it tends to like (1 - z/p)^ 1 )/* The estimate for G^ n) follows. 
Induction step for coalescing small branches u?,-: Assume that n > 1 and the assertions 
are true for all pairs (r, s), where r < n — 1. According to the representation (4.2), we have 



$(z,u)F n > u (u) 



RHS n (z, u)-J2 RHSn(z, ui{z)) 



det M{u; I) j-r u 



i=i 



detM 



n 

i=i 



U i 



(5.15) 



Ui 



The estimate for F n,t (wj) is computed by induction on t by taking t + 1 derivatives of this 
equation and evaluating at Wj. 



(t + l)F^( Wj (z)) = 



E 

fc=2 



t+1 

fc 

1 



ra.i+l-A; 



^'(z^z)) 



c 



(^•(z)) + 



— —-(t+i) 

i?F5 n (z,u;,-(z)) 



_9_ 
9« 



^'(z^z)) 
t+1 det.M(u;Z) 



detX 



n 

3=1 



U — Ui 



Ul 



(5.16) 



Observe that for z — > p and by (5.4) we have 



0_ 
du 



n 



u 



Ul 



U i 



O ((1 - z/p)- r /P) , u, G . . . , w p }, r < p - 1, 
O ((1 - z/p)(P- i yp) , ^ G . . . , w p }, r > p, 
O ((1 - z/p)(^)/P) ^ {«;!,... , r < p - 
O (1) , ui {wi, w p }, r>p. 
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Furthermore by the expansions (5.5) we see that 



¥ k \z,Wj) _ fo((l- z/pY^-^lv) , k<p-l 
&{z,Wj{z)) ~ \ O ((1 - z/p)-^-^/P) , otherwise. 



-~(t+i) 
By induction hypothesis, RHS n 



estimates (5.14) for RHS n (ui), i 

induction on t, by keeping track of the singular orders at z 



(z,Wj) = O ((1 - z / p )-(3n/2-2+(t+i)/p)) _ With the above 

■■ 1, . . . , c, we can prove the assertions for F n,t {wj{z)) by 

p. ' □ 



6 Limit distribut ion of excursion area. 



Before we apply the above result to prove the limit laws via coefficient asymptotics, we recall 
the Transfer Theorem 



Theorem 6.1 (Theorem VI. 3 in [E]). Let a G C \ {0, —1, —2, . . .} and F(z) be analytic in 
the open indented disc 

D(p, o-, </>) = {z G C | \z\ < a, z ^ p, | arg(z - p)\ < 4>} 



where < p < a and < (p < vr/2. // in i/te intersection of a small neighbourhood of p with 
D(p,a,4>) F(z) satisfies the condition 

F(z)~(l-z/ P r a (*— 



. n 



a-l 



i/ien [z n ] F(z) ~ p — r-, for n — > oo. 
Via) 



□ 



With Y(z,u) = W(z,u)Q(z,u) we have by (4.3) 



(0), _ det g k {Y) 
&k {Z) ~ detM ■ 



By lemmas 



5.2 



and 



5.3 



G. (z) has an expansion in non-negative powers of y/l — z/p. In order 
to access the coefficient of -y/l — z/p we regard formally as a function of the variable 
u\ into which we substitute the expansion u±(z) = r — fi\f\ — z/p + 0(1 — z/p). By Tailor's 
theorem we hence have 

G f\z) = Gf\p) - e^zTp + 0(1 - z/p) 



where 



det M mi det ~ det g k & det M 



e = lim ■ 

z— >p 

= lim 



detM' 



det.M M ^detM(z,u;y(z,u))\ u= 



2a ,oQ(/0, t) lim 



detM 2 
det M k A 



z^p det M 
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Here we used the relation (4.6) and the definition (5.8) of aoo- Furthermore, recalling that 
zQ(z,ui(z))S(z,ui(z) = Q(z,u\(z)), we have by Lemma 5.3 that 



RHS n (z,u 1 (z))=Q(z,u 1 (z))nF n - 1 ' 1 (u 1 (z)) + 0( K (l-z/p)- ( - 3n / 2 - 1 ^ (z — > p), (6.1) {eq: 



RHSestimate 



which dominates RHS n (ui(z)), i > 2. Therefore, with (4.3) and the leading order term given 
in 



5.3 we have for z 
1 



P 



ill 



lQ(z,u 1 (z))nF-^(n 1 (z)) d £^i + O ((1 - z/p)"^-!) 
1 Qn-1,1 n( s y det A^fc.i , 

(l_^/ p )3n/2-l/2( n _l)!^^' T ^™ det.M ^ 



(6.2) {eq:singGOn} 



For the random variables X m defined as in (2.1 ) we hence have by the Transfer Theorem the 
following asymptotics for the factorial moments, 



1 [z m ]G^\z) 



F(-l/2) 



»! [z™] G { °\z) 2a ,o(n - 1)! T(3n/2 - 1/2) 



in 



3n/2 



(6.3) {eq:asExMo} 



This shows that the ordinary and factorial moments are asymptotically equal and hence 
the moment convergence for the sequence ^m _3 / 2 X m . By the recursion equation (5.7), the 

normalised sequence C n ± 
hence have the moment convergence 



— n! fi "' is precisely the one from part i) of Theorem 



2.2 



We 



lim — E 

m— >oa 77,! 



-r(-i/2) 

-l 2 U n _n, 



.v^m 3 / 2 ; y r(3n/2-l/2) 
which implies the asserted weak convergence to Brownian excursion area [7j, once Theorem 



2.2 is established. This is done in Section 10 Our proof also shows that the same limit 



distribution arises, if we choose another Gk for the definition of X m , i.e. the distribution of 
area under meanders ending at a fixed altitude < c is asymptotically BE A. □ 



7 Limit distribution of Z m : negative drift 

We now consider the case when &(z, 1) 7^ for \z\ < p, which in the frame work of discrete 
meanders simply means that the drift 7 = S'(l)/S(l) is negative (and hence r > 1). It turns 
out, that the dominant role is played by the quantities F n,t (ui). 

Lemma 7.1. We have 

F^\l) = o{F^( Ul {z))) (z^p). 

Proof. Once the assertion is true for (n, 0), we can prove the statement for (n, t), t > 1 by 
induction on t as follows. Assume the statement to be true for all k < t and consider the tth 



derivative of (3.6) w.r.t. u evaluated at u = 1. By induction hypothesis, every (known) term 
is dominated by F n,t ~ 1 (u\{z)) for z — > p. Since $(z, 1) 7^ for \z\ < p, solving for F n,t (l) 
shows that it is also at most of order as F n,t ~ 1 (ui(z)) which is in turn o (F n, *(ui(z))) by 
Lemma 15.31 
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The assertion for (n, t) = (0,0) follows from the explicit representation (4.2) of F 0,0 (u 



and Lemma |5.3[ Now assume that the assertion is true for all (k,t) with k < n — 1. By 
the representation of F n,0 (u) obtained applying (4.2) with y(z,u) = RHS n (z,u), we see that 
F n -°(1) is O (F n - l ^( Ul (z))) as z — > p, which in turn is o(F n > ( Ul (z))) . □ 

Let again Y(z, u) = Q(z, u)W(z, u). We see easily that the singular behaviour of F ' (z, 1) = 
det Af(z, 1; Y(z, u))/(&(z, 1) det Af) is of square root type. Similarly to Section [6j we access 
the coefficient of \J\ — z/p by regarding F°'°(z, 1) formally as a function of the variable u\ 
into which we substitute the expansion u\ 

(z) = T - p^J\ - z/ p + 0(1 - z/p). By Tailor's 

theorem we hence have 



F°'°(z, 1) ~ F°>V l)\ ui=r - vPVl^fp (z — ► p), 



where r/ is given by 







lim — 



lim 

2— >/9 

lim 

Z— >p 



detAA(z,l ; y) gf^detAl 



det A* 2 £ r detJV(*,l;y) 



l)detA( ; 



detA1(l;l) ^ detA/(z, u; F) | 
l)detAf 2 



where we used the relation (4.7). Finally, we plug in definition (5.8) and get 

2a ,oQ(p,T) 



V 



lim 

/3$(p,l) 



det M(l; 1) 
det Af 



For the leading singular behaviour of F n ' (1) we recall (6.1) and see by virtue of Lemma 

IS n (z,ui(z)] 

nQ(z, ui(z)) 



7.1 



that RHS n (z,l) = o(RHS n (z,ui(z))). With the representation (4.3) 
F n -°(1) 



n! 



-1 

-l 



??! 



Q(p,r) 



det M(l;l) 
detAl 4 
detA((l;l) 



lim 



(n - 1)! (1 - z/p) 3 ™/2-i/2 l) det Al 



O ((1 - z/p)-( 3 «/2-l)) 



The moments of Z m are seen to be asymptoticaly equal to those in (6.3), namely 

1 [z m ]F n >°(z,l) 



E((Z m )J 



r(-i/2) 



[z m \ F°>°(z, 1) 2a ,o(n - 1)! T(3n/2 - 1/2) 



m 



3n/2 



The rest of the argument is as in the previous section. 



□ 



8 Limit distribution of Z m : zero drift 



As mentioned in the remark following Theorem 2.1 the convergence of the joint distribution 
of meander area and endpoint is a simple consequence of the weak convergence result in 
[17J. The formula (8.6) for the joint moments given below yields the recursion (2.14) in 
Theorem |2.2| If in turn we consider Theorem 2.2 as given, then the weak convergence to 
(BA4A, B me {l)) follows from the moment convergence (8.6), which is as a consequence of the 
functional equation. This is case 2 in Theorem 5.1 r = 1 and <&(z, 1) ~ Q(p, 1)(1 — z/p) as 
z — > p. 
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Lemma 8.1. The mixed derivatives evaluated atq = u = l have the following leading singular 
behaviour as z — > p. 



F n ' (z, 1) 



+ o 



1 



(1 - z / p fn/2+t/2+l/2 Z /p) 3n / 2+t / 2 

where the numbers b n ^ for (n, t) = (0, 0) and (n, t) = (0, 1) are given by 

bo,o = 2a 0i o and 6 ,i = ^ & o,o, 

for (n, t) with t > 1 by the recursion relation 

t(t - 1 



(* — > P), 



B 2 



b n ,t~2 + nb n -i,t+i, 



i.l) {eq:singFntl} 



1.2) 



.3) {eq:recbnt} 



and for (n, t) with t = and n > 1 



b n ,o = nb n _n + na n -ii. 



The numbers a n t are defined in Lemma 5.3, 



Proof. The calculations are again similar to those in the proof of Lemma 5.3 Let Y(z,u) = 
Q(z,u)W(z,u). We first give the proofs for (n,t) = (0,0) and (0,1). As for &o,0) we argue 
similarly as in the previous section that 



detA/"(z,l;Y) 
detM 



with 



p = lim 



Pp^fl-z/p + 0(l-z/p) (*—►/>), 
d detM(z,l;Y) 



>p du\ det.M 
Since &(z, 1) ~ Q(z, 1)(1 — zp) we have 

detAT(z, 1;Y) p 



F°>°(1) 



p (3 



Q(p,i) ^Y^zTp 

The formula for 6o,o follows from the simple observation that 



d 



detM(z,l;Y) 



d_ 
du 



det N(z,u;Y) 



U=U\ = 1 



For bo i we consider 



F o,i (1) 



& det Af(z, u; Y)\ u=1 detM(z, 1; Y)&{z, 1) 



&{z, l)detM 



$(z,l) 2 detA4 



As z — > p, the second summand is of order 0(1/ y/1 — zp) since det J\f(z, 1; Y) is of order 
0(y/l — z/ p) and $>'(z, 1) of order 0(1 — z/p). The first summand is precisely of order 1/(1 — 
z/p) and hence 



F U)1 (1) 



d_ 



detM(z,u;Y)\ u=ui=1 i 



90,0 



Q(p,l)detM 



1 — zp (3 1 — zp 



(z — > p). 
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The remaining cases are shown by induction. For the t > cases, the tth derivative of 
equation (3.6) can be rewritten as 

(1 - zS(-u))F n '*(u) =tzS'(z,u)F n ' t - 1 (u) 

Ct\ , n i t ,i (8.5) {eq:Fntl} 

+ 1 J zS"(z, u)F n ^- 2 {u) + zS{z, u)nF n -^ t+ \u) + R^u), 

where we collected the F k,l (u) with k + l<n + t or k + l<n + t and k < n — 2 and the terms 
involving \ j <n in R n ^{u). 

We now consider the cases t > 1. Upon setting u = 1, we see that 5'(z, 1) = 0(1 — z/p), 
and therefore the term zS'(z, l)i ?n,t_1 (l) is insignificant for leading order considerations. By 
induction, i ?n, *~ 2 (l) and F n ~ 1,t+1 (l) strictly dominate the terms in i? n; t(l) as z — > p, and 
the recursion relation for t > 1 follows by comparing leading singular terms, recalling that 
p = 5(1) and /3 2 = 25(1)/5"(1). 

In order to prove the relation for t = 0, we have to combine (4.2) with (3.6) again, leading 

to 



,m,0 



u) =^ V) ^nzQ(z, U )S(z,u)F n -^ l ( U ) - ^Q^^-bl^^lM 



where we used that Q(z,ui) = zQ(z,ui)S(z,ui). By induction hypothesis, for it = 1 and for 
z — > p only the first line of the previous equation is of interest. By Lemma 5.3 the terms 
with F n ~ 1,1 (ui) for I > 2 are also neglegible, which finally leads to 

™, 0m nzQ(z, 1)5(3, - n Q(z, u 1 )F n - 1 '\u 1 ) F^ 1 ' 1 (1) - F^ 1 ' 1 ( Ul ) 
i 1 1 ~ 777 77 rr ~ n- 



Q(z,l)(l-z5(z,l)) (l-z5(p,l)) 

The assertion follows from comparing leading singular terms. □ 

As above, an application of the Transfer Theorem gives the leading asymptotic behaviour 
of the joint factorial moments 

^■«*.). w,) = ^g^g ~ n uwnnl% + < 8 - 6 ' ~> 

This implies that factorial and ordinary moments are asymptotically equal. If n = 0, then 



E((fl- m ) t ) 



l^]i^(l) 6 ,mV 2 - 1 / 2 r(l/2) / V2^\ , R - 
~ 6o,om-V2r(t/2 + 1/2) = 1 — 1 2 1(1 " (8 - 7) 




where we used the easy fact that &o,t = bofitl/P 1 = &o,or(* + and the duplication formula 
for the Gamma function. The number 2*/ 2 r(l + t/2) is the tth moment of the Rayleigh 
distribution on [0, oo), given by the distribution function P([0, x)) = 1 — exp (— x 2 /2) . As the 
Rayleigh distribution is also uniquely determined by its moments [31], we have rederived the 
limiting distribution of the final altitude of a meander with zero drift as was previously done 
in pp. Notice that this was solely derived from the functional equation. 
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In Section 10 we shall infer from a separate treatment of the Bernoulli case that the 



bu,tP n+t 



uniquely determine the joint distribution of area and final point of 



numbers Q n ^ 

the Brownian meander. The Q n t are precisely those defined in Theorem 2.2 and hence the 
general "zero drift" case of Theorem 2.1 can be viewed as a corollary of Theorem 2.2 □ 



9 Concentrated Z m : positive drift 

Finally, we treat the case where $>(z, 1) does have a simple zero < zq < p. The following 
recursive description of the leading singular behaviour of F n,t (l) at zq can easily be proven 
along the same lines as in the previous section. The dominant singularity in this case turns 
out to be a pole, more precisely 



F n -'(1) 



e n ,t 



(1- z/zo) 2 ^ 1 ' 
where the numbers e n j are given by the recursion 

e n ,t = tje n j-i + ne n _i )t+1 . 



The constant 7 is given as in Theorem 5.1 



S'(z , 1) 



S(zo,l) + z f z S(z,l)\ z=zo 



and the initial conditions are 



eo,o 



detM{z,l;Y) 



Q(z,l)detM 



z=z 



and e n ^ = for n < or t < 0. In particular a direct computation gives e^o 
e 2,o = 6eo,o7 2 - An application of the Transfer theorem gives 



e ,o7 and 



E(Z„ 



ei,o«t 7 2 

— 777, 

eo, r(3) 2 



w ir 7 2\ e 2 ,o"i 4 7 2 4 w n2 



and hence the variance Y(Z m ) is o ( E (Z„ 
centration property follows. 



By Chebyshev's inequality, the claimed con- 

□ 



10 The Bernoulli case and proof of Theorem 2.2 



The Bernoulli case plays a significant role in our derivations. In particular, we have weak 
convergence of the Bernoulli excursions to Brownian excursion, and by a theorem of Drmota 
[9] moment convergence for any polynomially bounded functional (like area and endpoint), 
which we lack a proof for in more general cases. Furthermore, we want to exploit their 



combinatorial simplicity in the next section to prove Theorem 2.3 



The numbers a n j and b n t from Lemmas 5.3 and 8.1 are related to those from Theorem 

2 3n+tpn+t Gn t 



2.2 as follows, 



C n ,, 



n\t\ 



oo,o 



Qn,t 



(3 n+t b n , t 
n\t\ 6 n' 
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Equation (5.7) implies the recursion relation 



C n ,t = C n ,t-i + {t + 2)C n _ M+2 , C ,o = 1, C n ,t = 0, if t < or n < 0, 



(10.1) 



and similarly, from equation (8.3) follows 



Qn,t 



Qn,t-2 + (t+ l)Q n -l,t+l, if t > 1, 

Q n _i i i + 2-8- n C„_i,i, if t = 0, 



(10.2) 



with the initial values Qo,o = Qo,i = 1 an d Q n ,t = if t < or n < 0. The numbers C U) t and 
Qn,t do not depend on the parameters of the walk. We show with the help of the convergence 
results for the Bernoulli walks that 



C n - 



i.i 



8 n K n and Q nfl = Q n , 



(10.3) {eq:cn-llanddnO 



with K n and Q n as in Definitions 2.1 and 2.2 This settles i) of Theorem 2.2 and ii) follows 



from a Theorem of Petersen [25J , since the sequence M n t is by construction a sequence of joint 
moments (as a limit of such) and the projections are uniquely determined by the moment 
sequences M n fi and Mqj. 

Xf^ n and meander area Z, 



The limiting distributions of the excursion area X2 m 
^Bem £ Qr ^ e symmetric Bernoulli step set S = {—1, 1} with weigths s_i 
to equal BE A and BMA, respectively [32l El] ■ More precisely 



si 



m — 

1 are known 



X 



2 m 



(2m) 3 / 2 



B£A and 



Z n 



3/2 



BMA. 



(10.4) {limitBern} 



Recall that there are no Bernoulli excursions with an odd number of steps. Our above 
derivation expresses the same limit laws in terms of numbers a n ^ = OnT n an d &n,t = b^f n , 



cf. equations (6.3) and (8.6). Equating the respective expressions together with the above 



normalisation of a^f n and b^f n prove assertion (10.3) and hence the general case. 

One detail has still to be taken care of: The step set S = {—1,1} is not aperiodic, i.e. 
our general assumption under which we derived the singular behaviour of the F n ' (iti) and 
the F n,t (l) is not fulfilled. On the one hand, aperiodicity implies the domination property of 
the small branches. On the other hand it implies the uniqueness of the dominant singularity 
p. The domination property is trivially fulfilled, since there is only one small branch. What 
demands a thorough review is the occurrence of two dominant singularities. 

The step polynomial is S(u) = u^ 1 + u and the only small branch of the kernel equation 
1 — zS(u) = is 

VT-2V1 + 2Z 



ui(z) 



1 



2z 

having two dominant singularities z = p = 1/2 and z = —1/2, ui(±l/2) = ±r = ±1. The 
computations of the leading singular behaviour of F n,t (u\) at z = p = 1/2 are still valid, 
leading to a recursion (5.7) for the a^f n , with the parameter /3 = /3 Bern = y/2. Repeating the 
calculation at z = —1/2 also gives the leading singular term of the expansion about —1/2, 
with the very same coefficient a^f n . This implies that G^\z) has two dominant singularities 
at z = ±1/2 with leading terms as derived in equation (|6.2[), namely 



G. 



in). 



1 



2a Be m (n - 1)! (1 T 2z) 3 "/2-i/2 



(z— >±l/2), 
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This has to be taken into account in the process of coefficient asymptotics. The effect of the 
additional singularity eventually cancels out, as is seen when computing the asymptotics of 
the factorial moments, 



1 



in 



E((*2m)J 



i [z 2m ] gra 



IV. 



y 2m] 



Bern 
z n-l,l 



T(-l/2) (2 2m + (-2) 2m ) (2m) 3n / 2 - 3 / 2 



2a^ n (n 



1)! T(3n/2 - 1/2) (2 2m + (-2) 2m ) (2m)~ 3 / 2 



Since factorial and ordinary moments asymptotically coincide, equating the so obtained ex- 



pression for lim r , 



,E(((2 
-Cn-i,i 



m) 



-3/2 X 



2m 



Bern 
~ a n-lA 



with the one from Definition 

K n 2- n / 2 . K„ 







leads to 



23n— 1 ^Bern^n 



2a^ rn (n- 1)! 



iff! 



^Bern^n : 



which settles the first part of assertion (10.3). 



We proceed with the Bernoulli meander. The computation of the leading singular term 
of F n '*(l) at z = 1/2 leads to a recursion with /3 = /3 Bern = y/2. F n >*(l) is of course 

also singular at z = —1/2, but of neglegible order. One can show by induction that i ?n, *(l) 
0(F n - 1 > 1 (ui)) as z — > -1/2. We hence have by dO), (|10.4[) and Definition 



2.2 



1 



lim — E 

m— >oo 7l! 

which finally yields 



m 



3/2 



LBern 
°n,0 



r(i/2) 



Q n r(l/2)2-"/ 2 



n!^ c mr(3/2n + 1/2) Q r(3n/2 + 1/2) ' 



(10.5) 



n,0 



i.Bern 
°n,0 



(8 



Bern^n 



I^Bcrn 
0,0 



^Bern^n ' 



This settles the second part of assertion (10.3) and hence the proof of Theorem 2.2 



Remark. The case of a more general periodic step set can be treated similarly. The branches 



-2wijc/p Te 2Kij 



u\ and v\ coalesce in the points (z, u) = (pe 
has to argue more carefully for the separation property. 



0, . . . ,p — 1. Moreover, one 



11 Example: Column convex polygons 

For many polygon models on the square lattice the area distribution in a uniform fixed- 
perimeter ensemble has been studied and limit laws have been derived with the help q- 
algebraic functional equations satisfied by the area and perimeter generating function [26] [28| 
30j . This approach did not work out for column convex polygons due to lack a of "symmetry 
in horizontal and vertical edges". Recall that a polygon is called column convex if every 
intersection with a vertical line is convex. We will prove that the limiting area distribution is 

be a. 

Denote by F(z, q, u) the generating function of column convex polygons, in which z marks 
the half-perimeter, q the area and u the height of the rightmost column. These polygons can 



be built up by adding one column at a time, leading to a functional equation (4.1) with [3] 



u 2 (l — z) 2 zu 2 (2z — zu—l) zu 

S{Z,U) = - r^- r^, ro(z,U) = — ry- -, Tx{Z,U) ~ 



u) 2 (l — zu) 21 ' (1 — u) 2 (l — zu) ' ' 1 — u 
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and W(z,u) = Furthermore, G (z,q) = F(z,q,l) and Gi(z,q) = J^F(z, q,u)\ u=1 . 

Denominators can be cleared with Q(z, u) = (1 — u) 2 (l — zu) 2 and we see that <3?(z,u) = 
(1 — zS(z,u))Q(z,u), ti = Ti ■ Q, and W ■ Q are polynomials. $ is a polynomial degree 4 in 
u and since $(0,it) = (1 - u) 2 , there are two branches with Puiseux expansions 1 ± y/z + . . . 
about 0. The "+" branch plays the role of u\(z) and the other that of u^z). The point (p, r) 
is found by computing the resultant of &(z, u) and $>'(z, u). Critical values for z turn out to be 
0, ±1, 3+2^/2 and p = 3-2>/2. The value r = 1 + ^2 is the unique double zero of 3>(p, u), and 
S"(p, t) > 0. Within < \z\ < p none of u\, and U2 solves Q(z, u) = 0, since 3>(z, 1) = only 
for z = 0, 1 and 3>(z, 1/z) = only for z = 1. Both m and U2 can be continued analytically 
along (0,p) with u\(z) > 1 and 1*2(2) < 1. The latter is true close to z = and remains true 
throughout (0, p] for continuity reasons, since the equation <&(z, 1) = has only the solutions 
z = 0, 1. It follows that u\ has a square root singularity at p while U2 is analytic throughout 
\z\ < p. Finally we look at det M. = to(z, u\)ti(z, U2) — to(z, U2)ti(z, ui). We have 

det M. = z 2 (ui — 112)111112(1 — zui)(l — zu2)ip(ui,U2) 

and it remains to argue why (p(u\, U2) = (z(ui + U2)+z(z — 2)u\U2 — 2z + l) remains nonzero in 
< \z\ < p. Assume that ip(ui,U2) = for some zq and solve this equation for u%. Substitute 
the so-obtained expression for u\ in terms of Zo, U2 into the equation $(zq, ui(zq)) = 0, which 
yields a rational expression in z$ and U2 whose numerator ^(20,^2) is a polynomial in zo,U2 
of degree 4 in 112 with integer coefficients. So (zq, ^2(^0)) is a common solution of both 
<I>(z,u) = and *&(z, u) = 0. A possible zq is hence to be sought among the zeroes of the 
resultant of <3? and which are z = 0, zq = 0.167. . . < p and some other values of z with 
\z\ > p. The common solution of $(zo, u) = and ^(zq, u) = turns out to be 2.116 . . . > 1 
and hence cannot coincide with U2(z). So det.M 7^ in the domain in question and the 



assumptions of Section 5.1 are satisfied 



12 Unconstrained walks: joint distribution of final altitude 
and signed areas 

In this section we apply the result on the joint distribution of meander area and final altitude 
to get this joint distribution in a number of other cases. We have to restrict ourselves to 
the symmetric Bernoulli case as we want to use a decomposition which is restricted to that 
case, namely cutting the path whenever it crosses the x-axis, i.e. changes the sign. This 
yields a decomposition into positive and negative excursions which fails for general step sets. 
However, the rescaled Bernoulli walk converges to Brownian motion and the results carry 
over to that continuous limit making the considerations worthwhile. We consider the joint 
distribution of the final altitude and absolute area, area of the positive part and area of the 
positive and negative parts. The computations are pretty repetitive in the three cases, the 
most annoying of which is certainly the joint distribution of final altitude and the two signed 
areas. We demonstrate this latter one in more detail and only indicate the remaining. In this 
way rederive a number of formulae for the various Brownian areas as collected in the survey 
|18j . The joint moments with final altitude turn out to be simple generalisations of these. 

A bridge is a walk with steps from S = {—1, 1} that starts and ends at 0. Denote by 
B(z, q-, q+) its generating function, where z marks the number of steps, q + the positive and 
q^ the negative area. A bridge has either an even or an odd number of sign changes. In 
the even case it can be decomposed into a non-empty initial excursion and a sequence of 
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pairs consisting of a non-empty positive and a non-empty negative excursion. In terms of 
generating functions this reads after some simplifications 

B{z,q + ,q-) =: B(q + ,q-) = — — — - — — — -, 12.1 

Go(z, q+) + G (z, q-) - G {z, q + )G {z, <?_) 

in particular B(l, 1) = G^\z) ^2 — G^^z^j . Similarly, we decompose unconstrained walks. 

A walk is either empty, or it can be decomposed in a (possibly empty) initial excursion, fol- 
lowed by a (possibly empty) sequence of pairs of non-empty excursions as above, followed 
by a non-empty meander. Some simplifications finally yield for the generating function 
W(q + ,q-,u) := W(z, q + , g_, u) the representation 

W(q+,q-,u) = —— r— — — — — r-^ - B(z,q+,q-). (12.2 

G (z, q+) + G (z, q_) - G (z, q + )G (z, qj) 

For convenience we recall the singular behaviour of Gq (z) and F n,t (l) at z = 1/2 in the 
Bernoulli case, namely 

2-2v / 2VT^, G^\z)~ Q 00^ (z^l/2). (12.3) {GOGn} 

and 



n ] t l 2 a-n-t)/2Q nt 
(1 - 2z) 3n /2+*/2+l/2 • 



/n ^ttSW - (12-4) {FntBern} 



12.1 Signed areas and final altitude 

As in the previous sections, we introduce shorthand notation for the mixed derivatives eval- 
uated at q+ = q- = u = 1, namely 

W k,i,t = W k,i,t( z ) = (JL\ (JL\ fl \ W(l,l,l), (12-5) 



dq+ J \ dq- J \du 

analogously B k ' 1 and R k > 1 for R(q + ,q_) = (Go(z,q + ) + Go(z, qJ) — Go(z,q + )Go(z,q^))~ 1 is 
defined. 

A twofold application of (a univariate version of) Faa di Bruno's formula [8] yields 



dq + ) \dqj R ^+^ E R(q + ,q„)\M+\v\+l 



(12.6) {faadibruno4} 



where the sum runs over pairs of partitions (ir, rj) of a k- resp. l-set and the product over the 
blocks of the respective partitions. By | • | we denote the cardinality of a set. The omitted 
trailing terms are of smaller asymptotic order as q+ = q~ = 1 and z — > 1/2. Observing 
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that R(l, 1) = (G { q\z) [2 - G^'(z))) and G^\l/2) = 2 this yields together with the 
asymptotic representations (12.3) 



(0), 



k.l 



2 -(fc+0/2 



:Dt 



4^/2(1 - 2z) 3fc / 2 + 3Z /2+i/2 W 
where the numbers can by ( |12.6[ ) be expressed as 



(12.7) {Rkl} 



D^ = ^2(\n\ + \ V \y.l[K lbl \b\\K lcl \c\ 



x k y l 



7T,?7 



ben 



1-L>i^ n + ^ 



(12.8) {Dklpm} 



The numbers -D^ ; are related to the joint moments of the signed areas of the Brownian Bridge, 
see [TBI eqs- (251) and (252)] and [24J. Indeed, we further find that B k > l (z) ~ 4 J R fc >'(z) as 
z — > 1/2. Applying the product rule and regarding only the dominant terms we obtain 

i=o W j=o Uy 



i=0 

2-(*+'+*)/2jfeint! 
(1 - 2z) 3fc / 2+3Z / 2+ 



(12.9) {Wklt} 



1 n k,l,V 



where by equations (12.7), (12.3) and (12.4) we have the following representation of the L, lt 



1 - i (—1)* ■ 



(12.10) 



i=0 



j=0 



the factor (—1)* resulting from 



d\ n ( a V 

dqj \du 



\ F(z,q,u~ l ) 



= (_i)* (F n '\l)+t( y t-l)F n ' t - 1 (l) + . ..). (12.11) {Fntumimisl} 

■u=g=l 

We hence have moment convergence for the rescaled random variables, 



E 



Zm \ ( Z m \ ( H 



m 



3/2 



nr 



3/2 



m 



1/2 



fc!f!t!2-( fc +'+*)/ 2 
T(3A;/2 + 3//2 + t/2 + 1) 



r± 



In terms of generating functions we have the relation 

V- T ± uit E n , f >o Qn,t {x^ + !/«(-«;)*) 



(12.12) 



which for w = specialises to |181 eq. (275)]. Furthermore we obtain the joint distribution 
of the area of the positive part and terminal altitude by simply specialising to I = 0, which, 
for t = 0, gives [181 ec L- (298)]. The numbers D^ k = D^ are referred to as in |18| . 
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12.2 Absolute area and final altitude 



This is the case q+ = q- = q, in which the above formulae for the bridge and walk generating 
functions simplify to 

G (z,q) . . F{z,q,u) + F{z,q,u- 1 ) 
B{z > q) = 2-G (z, q y W{q > u) = 2-G (z,q) 

By an application of Faa di Bruno's formula we find similar to the above derivation that 



1 (d_ 

2 \dq 



B(z,l) 



d 
Oq 



G {z,q) 



9=1 



The numbers are given by 
1 



D,, 1 = r ^2W| 7 r|!n^|6|W 

vr ben 



k\ 



2- k / 2 D k 
2^2(1 - 2z) 3fc / 2+1 /2 ' 



l-2E„>iV 



En>0 ^n- 



where ir runs through the partitions of a fc-set and b through the blocks of ir. The numbers 
-Dfc are related to the moments of the absolute area of the Brownian bridge, cf. [181 eqs. (135) 
and (287)] and |19 t [2^ ], [33] . Comparing the generating functions of the numbers and D^- 

from ( |12.8 ) we see D k = Yli=o^k-a^ i^S eq. (254)]. The product rule then gives for 
even orders 2t 



a 

dq 



d_ 
du 



2t 



k=0 

n!(2i)!2- n / 2 -* 



' 2^2(1 _2z) 3fe / 2 +V2 



(1 - 2z) 3n / 2+t+1 

with L n ^t = J2 T k=oQn-k,2tDk- We have the generating function identity 

T n 2t K Y; n , t>0 Qn,2tX n y 2t 
2^ L n,2tX V = ■ 
n,t>0 



By singularity analysis we see that the joint moments of the rescaled random variables X n 
of absolute area and H m of final altitude converge to a finite limit, namely 



E 



X„ 



m 



3/2 



H r 



m 



1/2 



2f N 



n 



\ t \2-nl2-t 



T(3n/2 + t + l 



-L 



n,2t- 



For odd t, the term F n,t (l) cancels out and the joint moments are 0(m -1 / 2 ) and hence the 
limiting moments are 0, in accordance with the obvious invariance of the distribution of the 
both variables under reflection of the Bernoulli walk. Hence we set L n ^k+i = for k > 0. 
For t = we retrieve |18| formula (290)]. 
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13 Conclusion 



We have proven a limit theorem for random variables defined in terms of coefficients of power 
series which in turn are determined by a "g-shift" of a linear equation with one catalytic 
variable. The main application thereof are the limit distributions for the area under discrete 
excursions and meanders with an (almost) arbitrary finite set of steps, providing some new 
results and re-deriving some classical ones on walks on the integers. As another combinatorial 
application of this result we prove an area limit law for column convex polygons, which has 
so far only been done numerically. 

As a bonus we computed the joint distribution of several Brownian areas and the height 
of the endpoint in terms of their moments. In particular Theorem 2.2 yields new recursion 
relations for the moments of BE A and BMA. 

For meanders with a positive drift we know from the theory of random walks that the 
limit law for the centered and normalised area random variable is Gaussian, and we believe 
that this also holds in the more general setting. Our method, however, would demand a 
considerable modification to prove this, since knowledge of only the dominant term of the 
singular expansions is no longer sufficient. 

According to our assumptions 5.1 we need a good knowledge of the solutions of the kernel 
equation in order to make Theorem 5.1 effective, which is similar for [27\ Theorem 1.5]. Our 
assumption 5 on the determinant in Section 5.1 may be difficult to check, as the computations 
for column convex polygons may suggest. 
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